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Abstract 



Two and three dimensional random Ising models with a Gaussian distribution of couplings with 
variance J and non-vanishing mean value Jo are studied using the zero-temperature domain-wall 
renormalization group (DWRG). The DWRG trajectories in the (Jo, J) plane after rescaling can be 
collapsed on two curves: one for Jo/J>r c and other for Jo/ J < r c . In the first case the DWRG flows 
are toward the ferromagnetic fixed point both in two and three dimensions while in the second case 
flows are towards a paramagnetic fixed point and spin-glass fixed point in two and three dimensions 
respectively. No evidence for an extra phase is found. 

In some range of concentration of magnetic impurities in a non-magnetic host one observes a competi- 
tion between spin-glass and ferromagnetic order [jy . The phenomenon can be described by an Ising model 
in which couplings are distributed randomly with some non-vanishing mean value J and width J. This 
issue has been addressed by Sherrington and Kirkpatrick [|| in the case of a soluble infinite-range model. 
For the case of insoluble short range model Migdal-Kadanoff real-space rescaling || , Q] and computer 
simulations have been employed. It is established ||,|| that there is no finite-temperature phase 

transition in a two-dimensional (2d) Ising spin-glass (random Ising model with Jo = 0). McMillan (7| in- 
vestigated the 2d random Ising model for the general case, Jo not necessarily zero, using finite-temperature 
domain- wall renormalization group (DWRG). He found at small temperatures only two phases: ferromag- 
netic and paramagnetic. However, some later studies ]12] ], (T^| have found a finite temperature transition 
to a "random antiphase state", which has similar properties to a spin glass. Therefore more study of the 
subject is desirable. To my knowledge no DWRG study of three dimensional (3d) random Ising model 
(with J =/= 0) have been done. In this paper the 2d and 3d random Ising models are investigated using 
the zero-temperature DWRG In particular no evidence for the "random antiphase state" is found. 

The system is the Edwars- Anderson model fbl of an Ising spin glass with Hamiltonian: 



where Si — ±1, Jy 7^ only when lattice sites i,j are nearest neighbors. Each coupling constant Jy is 
an independent Gaussian distributed random variable with P(Jy) = exp(— ( Jy — Jo) 2 /2J 2 ). The 

lattices studied are d-dimensional square (d = 2) and simple cubic (d = 3) lattices of linear size L. The 
zero-temperature DWRG method ^]-@ was used by computing the ground-state energies for periodic 
and antiperiodic boundary conditions (BC) in one direction, with free BC in the other d — 1 directions. 
The difference AE(L) = E p (L) — E ap (L) is the domain- wall energy which may be interpreted as an 
effective coupling constant on scale L. AE(L) has sample to sample fluctuations and one can define the 
mean Jq(L) and the width J{L) at scale L by 




(1) 



J (L) =< AE(L) > 



J(L) = \/< (AE(L)- J (i)) 2 >, 



(2) 



where < . . . > is the average over samples with different realizations of the coupling constants { Jy}. The 
ground state energy is found by simulated annealing fl5|| . Each annealing was start from random spin 
configuration or T = oo. The temperature was reduced in n s steps from Th — 5 to T) — 0.05 by a factor 
10~ 2 /" s . At each intermediate temperature T n = T^IO -2 ™/™ 3 the spins are updated by a single trial flip 
using the Metropolis algorithm. As a single annealing is not guaranteed to reach the minimum energy 
state, this procedure was repeated n a times, with the required CPU time proportional to n a x n s . It is 
found empirically that the lowest energies for fixed CPU time are found when n a ~ n s . This annealing 
schedule is substantially better than a large number of instantaneous quenches, for the same CPU time. 
To take maximum advantage of vectorization 64 annealings were performed simultaneously. To estimate 
the necessary values of n a and n s a few samples are tested until increasing n a and n s no longer leads 
to lower energy. Even so there is no guarantee that theses values of n a and n s are sufficient for all 
of the samples. To estimate the errors due to this, two runs for the same set of one hundred samples 
with periodic BC but with different sequences of random numbers for the Metropolis spin updating 
are performed. The two energy minima Ei,E2 found in the two nominally identical annealings were 
recorded. If true absolute minima are found than E\ — E2 but, if the true minima are not reached, then 
< (Ei — E2) 2 > is a measure of the error due to not finding the exact minimum. The numbers n a and 
n s are then increased until this error is less than statistical error due to the finite number of samples, i.e. 
until SE/E < A -1 / 2 . All data is obtained for N — I0 4 and the number of annealings n a and number of 
steps n s ranged from n a — 64, n s = 10 for the smallest sizes (L=2) to n a = 192, n s = 200 for the largest 
(L = 9 in d = 2, L = 5 in d = 3). 

The results of the simulations arc shown in Fig. 1 as renormalization group flows in the (Jo, J) plane. 
All initial couplings Jo(l) and J(l) are on the Jo + J = 1 line from which the trajectories, indicated by 
dashed lines, start. Next symbols on the trajectories correspond to L = 2, 3 

All the trajectories collapse into two curves by L - independent rescaling 

J (L) -> A(r)J (L) 

J(L) - \{r)J(L). (3) 

In Eq.3 the rescaling factor A(r) depends only on the ratio r = Jo/ J, and upon performing this 
rescaling all the flows of Fig.l collapse on to the flows of Fig. 2. That the rescaling of Eq. 3 works means 
that (at zero temperature) on any length scale L the system is described by the Hamiltonian of Eq. 1 
and no other terms are generated by the RG transformation. 

The flows for r > r c seem to be flowing to a ferromagnetic fixed point with J(L)/ Jq(L) decreasing and 
Jo(L) increasing with increasing L and one may speculate that the L = 00 fixed point is at J(oo)/ Jo(oo) = 
0,Jo(oo) = 00. For r < r c , on the other hand, the flows would seem to be different in d — 2 and d = 3. 
In d = 2 the flow is shown in Fig. 2(a) and it seems that as L increases both Jq(L) and J(L) ultimately 
decrease to zero with J(L)/ J (L) — ► 00. This would be interpreted as a paramagnetic fixed point for 
any finite temperature and a spin-glass fixed point at T = 0. In d = 3 (see Fig. 2(b)) the couplings 
Jo(L) — > and J(L) — * 00 which has the interpretation of a pure spin-glass fixed point. This spin-glass 
at T — is believed to survive at low temperatures T < T c in d — 3 fiofl . Note that the critical value 
of r = Jo/ J is r 2d = 1.05 and r^ d = 0.7 which is in fairly good agreement with the values r\ d = 1.2 and 
r^ d = 0.65 obtained by a Migdal-Kadanoff RG method The value of r 2d is in agreement with the 
result of McMillan [0, who extrapolated his finite-temperature DWRG data to T = and got r\ d = 1.04. 

In a study of a related two-dimensional Ising model jL2j, jl3) with a bimodal distribution P(Jy) = 
p5(Jij — 1) + (1 — p)S(Jij + 1) it was suggested that Jq(L) and J(L) scale as 

Jo (J) - L a 

J(L) ~ L d , (4) 

where for some range of p, a < and a > (sec p2[ ,Fig. 9). This was interpreted as evidence for a 
"random antiphase" state with properties similar to those of spin-glass. The raw data of Fig. 1(a) seems 
to support the conjecture of Eq.4 as there are trajectories along which Jq(L) and J{L) seem to scale in 



2 



this way. However, the simple additional rescaling of Eq.3 implies that the apparent a > is a finite size 
effect which disappears at large L. 

In conclusion this T = DWRG study of random Ising models with < Jy >= Jo > indicates that, 
at low temperatures, there are only paramagnetic and ferromagnetic phases in d = 2 and spin-glass and 
ferromagnetic in d = 3. No sign of any other phase is seen. Unfortunately, simulated annealing becomes 
ineffective for large system sizes L, and computer limitations restricted L to be quite small. However, 
recently it was found (Ref. jl6| ) that so-called branch and cut method allows one to study severall times 
larger systems than simulated annealing. Implementing this method for this particular problem may be 
a subject of future investigations. 
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Figure 1: Domain wall rcnormalization group trajectories for two (a) and three (b) dimensional random 
Ising models in the mean value (Jo) - variance (J) plane. The trajectories, indicated by dashed lines, 
start on the line Jo + J = 1. Direction of flows is indicated by arrows on some of the trajectoriies. The 
values of initial couplings Jo and J are indicated near corresponding symbols. Statistical uncertainty on 
each point is about the size of the symbol. 
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Figure 2: Rcscalcd, using Eq.3, renormalization group trajectories of d = 2 (a) and d = 3(b) random Ising 
models. Instead of Jo and J as in Fig.l their ratio is indicated near corresponding symbols. Trajectories 
for Jo/ J > r c {r1 d — 1.05, r zd = 0.7) flow to the ferromagnetic fixed point, while those for J / J < r c 
flow to paramagnetic fixed point in 2d and spin-glass fixed point in 3d. 
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